Time-domain models for power system stability and unbalance

It is an important and difficult challenge to protect modern interconnected power system from blackouts. Applying advanced power system protection techniques and increasing power system stability are ways to improve the reliability and security of power systems. Phasor-domain software packages such as Power System Simulator for Engineers (PSS/E) can be used to study large power systems but cannot be used for transient analysis. In order to observe both power system stability and transient behavior of the system during disturbances, modeling has to be done in the time-domain. This work focuses on modeling of power systems and various control systems in the Alternative Transients Program (ATP). ATP is a time-domain power system modeling software in which all the power system components can be modeled in detail. Models are implemented with attention to component representation and parameters. The synchronous machine model includes the saturation characteristics and control interface. Transient Analysis Control System is used to model the excitation control system, power system stabilizer and the turbine governor system of the synchronous machine. Several base cases of a single machine system are modeled and benchmarked against PSS/E. A two area system is modeled and inter-area and intra-area oscillations are observed. The two area system is reduced to a two machine system using reduced dynamic equivalencing. The original and the reduced systems are benchmarked against PSS/E. This work also includes the simulation of single-pole tripping using one of the base case models. Advantages of single-pole tripping and comparison of system behavior against three-pole tripping are studied. Results indicate that the built-in control system models in PSS/E can be effectively reproduced in ATP. The benchmarked models correctly simulate the power system dynamics. The successful implementation of a dynamically reduced system in ATP shows promise for studying a small sub-system of a large system without losing the dynamic behaviors. Other aspects such as relaying can be investigated using the benchmarked models. It is expected that this work will provide guidance in modeling different control systems for the synchronous machine and in representing dynamic equivalents of large power systems.

Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

Synthesis, structure, and stability of metal-organic frameworks

Metal-organic frameworks (MOFs) obtained much attention because of their unusual structures and properties as well as their potential applications. This dissertation research was focused on (1) the effects of synthesis conditions on the structures of MOFs, (2) the thermal stability of MOFs, (3) pressure-induced amorphization, and (4) the effect of high-valent ions on the structure of a MOF. This research demonstrated that the crystal structure of MOF-5 could be controlled by drying solvents. If the vacuum solvent is dimethylformamide (DMF), the crystal structure of MOF-5 is tetragonal. In contrast, if the DMF is displaced by CH2Cl2 before the vacuum, the obtained MOF-5 occupies a cubic structure. Furthermore, it was found that the tetragonal MOF-5 exhibited a mediate surface area (300-1000 m2/g). The surface area of tetragonal MOF-5 is also dependent on Zn(NO3)2/H2BDC (H2BDC: terephthalic acid) molar ratios used for its synthesis. The optimum ratio is 1.38, at which synthesized tetragonal MOF-5 exhibits the highest crystallinity and surface area (1297 m2/g). The thermal stability and decomposition of MOF-5 were systematically investigated. The thermal decomposition of cubic and tetragonal MOF-5s resulted in the same products: CO2, benzene, amorphous carbon, and crystal ZnO. The thermal decomposition is due to breaking carboxylic bridges between benzene rings and Zn4O clusters. Identifying structural relationships between crystalline and noncrystalline states is of fundamental interest in materials research. Currently, amorphization of solid materials at ambient temperature requires an ultra-high pressure (several GPa). However, this research demonstrated that MOF-5 and IRMOF-8 can be irreversibly amorphized at ambient temperature by employing a low compressing pressure of 3.5 MPa, which is 100 times lower than that required for amorphization of other solids. Furthermore, the pressure-induced amorphization (PIA) of MOFs is strongly dependent on the changeability of bond angles. If the geometric structure of a MOF can allow bond angles to be changed without breaking bonds, it can easily be amorphized by compression. This can explain why MOF-5 and IRMOF-8 can easily be amorphized via compression than Cu-BTC. It is generally recognized that zeolitic imidazolate frameworks (ZIFs) occupy much higher stability than other types of MOFs. The representative of ZIFs is Zn(2-methylimidazole)2 (ZIF-8) exhibiting high-decomposition temperature and high chemical resistance to various solvents. However, so far, it is still unknown whether the high stability of ZIF-8 can be challenged by ions, which is important for its modification by doping ions. In this research, we performed aqueous salt solution treatment on ZIF-8, and the results showed that anions (Cl¯ and NO3¯) in a solution exhibited no effect on the crystal structure of ZIF-8. However, the effect of cations (in a solution) on structure of ZIF-8 strongly depends on the cation valences. The univalent metal cations showed no effect on the structure of ZIF-8, whereas the bivalent or higher-valent metal cations caused the collapse of ZIF-8 crystal structure. Therefore, structure stability of ZIF-8 is considered when it is subjected to the application, in which high-valent metal cations are involved.

STUDY OF SPARK IGNITION ENGINE COMBUSTION MODEL FOR THE ANALYSIS OF CYCLIC VARIATION AND COMBUSTION STABILITY AT LEAN OPERATING CONDITIONS

A fundamental combustion model for spark-ignition engine is studied in this report. The model is implemented in SIMULINK to simulate engine outputs (mass fraction burn and in-cylinder pressure) under various engine operation conditions. The combustion model includes a turbulent propagation and eddy burning processes based on literature [1]. The turbulence propagation and eddy burning processes are simulated by zero-dimensional method and the flame is assumed as sphere. To predict pressure, temperature and other in-cylinder variables, a two-zone thermodynamic model is used. The predicted results of this model match well with the engine test data under various engine speeds, loads, spark ignition timings and air fuel mass ratios. The developed model is used to study cyclic variation and combustion stability at lean (or diluted) combustion conditions. Several variation sources are introduced into the combustion model to simulate engine performance observed in experimental data. The relations between combustion stability and the introduced variation amount are analyzed at various lean combustion levels.